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We show that inhnite Matrix Product States (MPS) constructed from conformal field theories 
can describe ground states of one-dimensional critical systems with open boundary conditions. To 
illustrate this, we consider a simple inhnite MPS for a spin-1/2 chain and derive an inhomogeneous 
open Haldane-Shastry model. For the spin-1/2 open Haldane-Shastry model, we derive an exact 
expression for the two-point spin correlation function. We also provide an SU(n) generalization of 
the open Haldane-Shastry model and determine its twisted Yangian generators responsible for the 
highly degenerate multiplets in the energy spectrum. 
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Introduction. — For a long time, it has been known 
that the main curse of quantum many-body theory is 
the exponential growth of the Hilbert space dimension 
with respect to the number of constituting particles. In 
the last decades, the study of entanglement has signih- 
cantly alleviated this curse, at least to some extent, by 
recognizing the fact that only a tiny corner of the Hilbert 
space, with small amount of entanglement, is pertinent 
for the low-energy sector of Hamiltonians with local in¬ 
teractions. This deep insight lies at the heart of tensor 
network states [1], a family of trial wave functions de¬ 
signed for efficiently representing the physically relevant 
states in the tiny corner. The best known instance among 
them is the Matrix Product States (MPS) in one spatial 
dimension, described in terms of local matrices with finite 
dimensions. Their entanglement entropies are bounded 
by the local matrix dimensions, which are nevertheless 
sufficient for accurately approximating gapped ground 
states of one-dimensional (ID) local Hamiltonians [1113]. 
This discovery not only provides a transparent theoreti¬ 
cal picture for real-space renormalization group methods 
but also leads to a recent complete classification of 
all possible ID gapped phases [6]-[8]. 

For ID critical systems, the low-energy physics is usu¬ 
ally described by conformal field theories (CFT). Their 
ground-state entanglement entropies exhibit unbounded 
logarithmic growth HE] with respect to the subsys¬ 
tem size, indicating the deficiency of a usual MPS de¬ 
scription. To overcome this difficulty, infinite MPS, 
whose local matrices are conformal fields living in an 
infinite-dimensional Hilbert space, have been introduced 
in Ref. [12]. The lattice sites for the infinite MPS lo¬ 
cate on a unit circle, embedded in a complex plane. This 
construction shares conceptual similarity to Moore and 
Read’s approach m of writing 2D trial fractional quan¬ 
tum Hall states in terms of conformal blocks. For a vari¬ 
ety of examples [T^ [T4H2Q] . the infinite MPS (as well as 
their parent Hamiltonians) have been shown to describe 
critical chains with periodic boundary conditions (PBC) 


and, furthermore, their critical behaviors are often re¬ 
lated to the CFT whose fields are used for constructing 
the wave functions [21]. In this sense, the infinite MPS 
introduced in Ref. m provide a systematic way of find¬ 
ing lattice discretizations of CFT. 

In this Rapid Communication, we show that the infi¬ 
nite MPS ansatz can describe ground states of ID criti¬ 
cal systems with open boundary conditions (OBC), thus 
complementing the PBC case in Ref. M- Unlike hulk 
CFT for periodic chains, open critical chains are instead 
described by boundary CFT. Taking a spin-1/2 chain as 
an example, we show how the infinite MPS with an image 
prescription allows us to derive an inhomogeneous open 
Haldane-Shastry model, including the original spin-1/2 
open Haldane-Shastry models [22] |23] as special cases. 
Within the new formalism, an exact expression for the 
two-point spin correlator of the spin-1/2 open Haldane- 
Shastry model is obtained. This, together with numerical 
results for the entanglement entropy, is in perfect agree¬ 
ment with the theoretical predictions based on boundary 
CFT, which thus confirms that our infinite MPS with the 
image prescription is suitable for describing open critical 
chains. The open infinite MPS construction is readily 
applicable to any boundary CFT for finding their lat¬ 
tice discretizations. As a further example, we derive an 
SU(n) generalization of the open Haldane-Shastry model. 
We characterize its full spectrum and also determine the 
twisted Yangian generators responsible for the highly de¬ 
generate multiplets in the energy spectrum. 

Infinite MPS and parent Hamiltonian. — Let us con¬ 
sider a spin-1/2 chain located on the upper unit circle 
in the complex plane, with L lattice sites and complex 
lattice coordinates Zj = {j = and Oj G [0, tt] 

Vj), see Fig. B a). We denote by 5'^ (a = 1,2,3) the spin- 
1/2 operators at site j. The local spin basis is defined by 
|sj), where Sj = ±1 (twice of the Sj projection value). 
For each site, we introduce its mirror image in the lower 
unit circle, e.g., site j has an image j, with complex co¬ 
ordinate Zj = zfi Following Ref. [12], the wave function 
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FIG. 1: (Color online) Schematic of an open chain in the 
upper complex plane. The lattice sites and their mirror im¬ 
ages locate on the upper and lower unit semicircles, respec¬ 
tively. They are symmetric with respect to the real axis. The 
two (brown) lines denote the chord distances \zi — Zj\ and 
\zi — Zj\, respectively, (a)-(c) denote the three uniform cases: 
(a) type-I: 0j = f (j- |); (b) type-II: 6>j = (c) type-III: 

0- — i 

~ 2L + lJ- 


is written as a chiral correlator of CFT fields: 


^-(51,... ,sl) = (1) 

where (uj) = Xj • :(:...: denotes normal 

ordering) and uj = [zj Zj)/2 (i.e., Uj is the coordinate 
of the “barycenter” of j and j on the real axis). Here, 
(l){u) is a chiral bosonic field from the c = 1 free boson 
CFT, and Xj = j even, respectively. 

Evaluating the chiral correlator in Q yields a Jastrow 
wave function 

4'(si, . . . , Sl) = <5^6*5 - It; ) ^ (2) 

j<l 

where (5^ = 1 if sj = 0 and zero otherwise (note 

that L must be even for ensuring a nonvanishing wave 
function). From the explicit form it is transparent 
that the sign factor (originated from Xj) is the “Marshall 
sign”, since the Jastrow product in ([^ is positive. 

As shown in Ref. m, the infinite MPS ([T]) with coor- 
dinate choice Uj = Zj = i.e., the case of equidis- 

tantly distributed spins on the whole unit circle, yields 
the ground state of the SU(2) Haldane-Shastry model 
[24l [25], which is a paradigmatic spin-1/2 chain with 
PBC. 

Now we demonstrate that our infinite MPS 0 with the 
image prescription, Uj = {zj^z^)/2^ describes a spin-1/2 
chain with OBC. Let us first derive a parent Hamilto¬ 
nian for which 0 is the exact ground state. Based on 
the CFT null field techniques, it was shown [14] that the 
decoupling equations satisfied by 0 lead to a set of oper¬ 
ators annihilating the wave function ([^, = 0 Vi, a, 

where Cf = | + iSabcSfS]) and Sabc is 

the Levi-Civita symbol [we assume summation over re¬ 
peated indices and use the convention that 
sum over j, whereas is the sum over both i and j]. 

When adapting to our present OBC setup, we consider 
the operators A“ = | +Wij){S^ + iSabcSiS^), 


where Wij = {zi + Zj)/{zi — Zj) and which also anni¬ 
hilate the wave function |4^), since = {zi — z^)Cf. 
The parent Hamiltonian for ([^ is then defined as = 
I f + E, where Si ■ Sj is the 

total spin operator and E = ^ + '^ij) “ 

After some algebra [26] , we arrive at a long-range Heisen¬ 
berg model 


H = 


E 


'y-. I 2 


Wj-jiCj+Cj) 

12 


{Si ■ Sj) 


WjjjCj - Cj) 

12 

( 3 ) 


with ground-state energy E, where Cj = wjj + 

+mj)- 

Three choices of the lattice coordinates deserve special 
attention (see Fig. (i) type-I: Oj = f (j — (ii) 
type-II: 0j = j^j; ^i) type-III: 0j = For these 

three cases (termed as uniform cases afterwards), one ob¬ 
tains Wij{ci—Cj)^w^j{ci^Cj) = 0, 4, and 2, respectively. 
Accordingly, the parent Hamiltonians, after removing the 
(unimportant) total spin operator and constant terms 
in 0, have purely inverse-square exchange interactions 
(between the spins and also their images), which coincide 
with the open Haldane-Shastry models first introduced in 
Refs. [22l[23]. These uniform models are integrable and 
have highly degenerate multiplets in their energy spec¬ 
trum [22] [23] , similar to their periodic counterpart [28] , 
see Fig. |^for the full spectrum of the open and periodic 
Haldane-Shastry models with L = 6. We postpone the 
discussion of this degeneracy until presenting the SU(n) 
generalization of these models, where a unified treatment 
is possible. The Hamiltonian 0 with lattice coordinates 
other than the three uniform cases is an inhomogeneous 
generalization of the open Haldane-Shastry models and 
does not exhibit the huge degeneracy in the spectrum. 

Spin eorrelator .— A nontrivial application of the in¬ 
finite MPS formulation is that, for the wave function 
0, the spin correlation functions can be computed eas¬ 
ily. Since = 0, one has (^| = 0 

and ~ ^ which lead to a set 

of linear equations relating two-point correlators Cij + 
= -f 13, where Cij = ■ 

5'j|4^)/(4^|4^). These equations are sufficient for comput¬ 
ing the two-point spin correlators for arbitrary choices of 
Oj (both inhomogeneous and uniform cases). The gen¬ 
eralization to arbitrary higher-order spin correlators is 
rather straightforward. 

Most remarkably, for the type-I uniform case, these 
linear equations allow us to find an analytical expression 
for the two-point spin correlator [26] 


L/2 p-1 


Cij = 




3(—1)* ^ sinOisinO. 
L{cos0i-cos0j) 

X cos 2q0j — cos 2q0i cos(2p — 1)0j 


( 4 ) 
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FIG. 2: (Color online) The energy spectrum of the three types 
of spin-1/2 open Haldane-Shastry models and the spin-1/2 pe¬ 
riodic Haldane-Shastry model {H = '^i^j sin2[f(z-j)/L] ) 

L = 6. All four models have highly degenerate multiplets 
in their energy spectrum. While the first excited states of 
the periodic model are degenerate singlet and triplet (due to 
two free spin-1/2 spinous), the open models do not have this 
degeneracy, indicating the importance of the boundary effect. 




FIG. 3: (Color online) Two-point spin correlators of the wave 
function in the type-I uniform case with L — 100. The 
blue circles are the exact results from Q , and the red crosses 
are fits with theoretical predictions based on the SU(2)i 
WZW model with free boundary condition (see text), (a) 
Two spins at lattice sites 50 and 50 + r are far from the 
boundary, (b) One of the spins lives at the boundary (the 
first spin). For (a) and (b), the first four points are excluded 
when computing the fits, since the theoretical predictions are 
valid for large r. (c) Two spins are nearest neighbors. 

with 

f 1 p=l,5 = o 

gp, = { p>hq = o . (5) 

p>i,9>o 

In Fig. various correlators from @ are compared 
with the theoretical predictions [29] based on the SU(2)i 
Wess-Zumino-Witten (WZW) model with free bound¬ 
ary condition. When two spins at sites j and j + 


r are both far from the boundary, one expects that 
the correlator Cjj-^r recovers the result for PBC [30] . 
Cjj+r ~ ai(-lG[vsin(f£)] oc (-l)7r for large r, 
where ai is a constant. However, if one of the two 
spins (say, the one at site j) is very close to the bound¬ 
ary, the theory developed in Ref. [29] predicts Cj j-^r — 
a 2 (-l)’-[^sin(y)]V 2 /[^sin(f £)]2 oc (a^: 

nonuniversal constant) with a boundary critical expo¬ 
nent ?7 = 3/2 that differs from 77 = 1 in the bulk. For 
the correlator between nearest neighbors, it was pre¬ 
dicted [Sll [ 32 ] that = as + a 4 /[^ sin(^)]^ + 

a 5 ( —1)-^/[^ sin(^)]^, where K is the Luttinger param¬ 
eter, K = 1/2, and as, a 4 , as are constants. We treat the 
nonuniversal constants ai,...,a 5 as fitting parameters 
and find excellent agreement between the exact result 
Q and the SU(2)i WZW predictions (see Fig. [^. 

Entanglement entropy .— To provide further support 
that the wave function (§ is relevant for open criti¬ 
cal chains, we numerically compute the Renyi entropy 
5 '( 2 )(/) = — InTrp^^ via Monte Carlo method [121 133 5 
where pi is the reduced density matrix of the first I spins. 
In Fig. 1^ we plot for the wave function ([^ in the 

type-I uniform case with L = 100. For open spin-1/2 
chains described by the SU(2)i WZW model with free 
boundary condition, one expects the Renyi entropy to be 

m 


5(2) (0 = ^ In 



+ C2 + 


(- 1)72 

[fsin(y)]f 


(6) 


with central charge c = 1, Luttinger parameter K = 1/2, 
and C 2,/2 nonuniversal constants. Fixing c = 1 and 
K = 1/2 and treating C 2 , /2 as fitting parameters, the nu¬ 
merical results are in good agreement with the theoretical 
prediction (see Fig. |^. For the type-II and type-III uni¬ 
form cases, we have verified via Monte Carlo simulations 
that their Renyi entropies also agree with suggesting 
that they all belong to the SU( 2 )i WZW model with free 
boundary condition. 

SU(n) generalization. — As a further application 
we generalize the above SU( 2 ) example to the SU(n) 
case. For the SU(n)i WZW model, the infinite MPS 
have been proposed in Refs. [m nH]. Here we take 
in all sites SU(n) spins transforming under funda¬ 
mental representations, with local basis denoted by 
\a) {a = l,...,n). Following Ref. [17], the CFT 
fields for defining the infinite MPS 0 are given by 
A^{u) = Ka : where rfio^ is a (n — 1 )- 

component vector denoting the fundamental weight 
of \a) (e.g., mi ,2 = (±l,l/\/3) and m 3 = (0,2/\/3) 
for SU(3), see [II])? is a vector of n — 1 chiral 

bosonic fields, and hia is a Klein factor, commuting 
with vertex operators and satisfying 
Evaluating the CFT correlator 0 , the SU(n) wave func¬ 
tion takes a simple Jastrow form, ^su(^)(ai ,... ^ai,) = 

sgn(k^\ • • • ,4%’ • • 
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monodromy matrix found for the spin- 1/2 open Haldane- 
Shastry models [23] to the SU(n) case. Through a third- 
order expansion of the monodromy matrix [26| , we obtain 
the nontrivial conserved charge responsible for the SU(n) 
open Haldane-Shastry models 0 

--i2 X] [Wjk+W^i,) 

~ '^ij)tk^jkPij ^ ( 8 ) 


FIG. 4: (Color online) Renyi entropy of the wave func¬ 

tion 0 in type-I uniform case with L = 100 as a function of 
the subsystem size 1. The blue circles (with errorbars) are 
obtained from Monte Carlo simulations and the red crosses 
are fits based on the theoretical prediction ® of the SlR2)i 
WZW model. The fit is computed with I G^0,90], as ^ is 
valid for large subsystem sizes. 


Uj^cxiocj (gg] 2 : signature of a permutation), where 
{k = 1 ,..., L/n), for a given configuration |ai,..., 
is the position of the kth spin in the state \a). 

Following a procedure similar to the SU( 2 ) case [26|, we 
obtain a two-body parent Hamiltonian for Tsu(n )5 H = 
^ r 1 I 1 Wij{Ci-Cj)+Wi-j{Ci+Cjy 






4(n+l) 


{ti ' tj), 

where {a = l,...,n^ — 1 ) are SU(n) generators in 
the fundamental representation, normalized as tr(t^t^) = 
^6ab- The three uniform choices of Oj, very much the 
same as the SU( 2 ) cases, bring the parent Hamiltonian 
into SU(n) open Haldane-Shastry models 




SU(n) 


= E 


1 


1 




{ti-tj) (7) 


with purely inverse-square interactions. 

Motivated by the SU( 2 ) result [23], we have numer¬ 
ically observed that the full spectrum of the SU(n) 
open Haldane-Shastry model 0 is described by the for¬ 
mula, i^su(n)|{^i}) = [^0 ^T.fLiE{mi)]\{mi}), where 
Eo = - mi - 1)] and E{m,) = 

l{m‘f - lN‘^) {N = 2L, 2L y 2, and 2L + 1 for the 
three uniform cases, respectively), M is an integer sat¬ 
isfying M G [0, and rrii are distinct integer/half¬ 

integer rapidities (1 < rrii <1/ — l, 2 <m^ <1/, and 
\ < rrii < L — ^ ioi each individual uniform case), sat¬ 
isfying the generalized Pauli principle which is the same 
as that for the SU(n) Haldane-Shastry model with PBC 
[3l|35]: only those sets {mi,..., tyim} without n or more 
consecutive integers/half-integers are allowed [28] . 

Twisted Yangian. — Our numerical results also indi¬ 
cate that the “supermultiplet” structure in the spectrum, 
which already shows up in the SU( 2 ) case (see Fig. 2 ), 
persists in the SU(n) open Haldane-Shastry models (7). 
To explain this degeneracy, we slightly generalize the 


where w^o = {^k + l)/(^/c ~ 1 ): Eij swaps the spin states 
at site i and j (more explicitly, Pij = 2ti -tj + ^) and 
7 i and 72 are given by (i) type-I: 71 = 0, 72 = 5 ; (ii) 
type-II: 71 = 0, 72 = (iii) type III: 71 = 1, 72 = |, 
respectively. The conserved charge and the total spin 
both commute with Q , but does not 
commute with the SU(n) Casimir operator 
This explains the appearance of degenerate eigenstates 
with different SU(n) representations. As the monodromy 
matrix relevant for these models (with open boundaries) 
satisfies the reflection equation [36], the algebraic struc¬ 
ture of the SU(n) open Haldane-Shastry models 0 is 
the twisted Yangian m- Thus, the conserved charges 
and T" form the lowest twisted Yangian generators. 

Conelusions. — In this Rapid Communication, we have 
shown that infinite MPS with the image prescription 
are relevant for ID critical chains with OBC, by pre¬ 
senting a spin- 1/2 example, as well as its SU(n) gen¬ 
eralization. We have constructed inhomogeneous open 
Haldane-Shastry models as their parent Hamiltonians, 
including the three open Haldane-Shastry models as spe¬ 
cial uniform cases. For the type-I spin- 1/2 open Haldane- 
Shastry model, an exact expression for the two-point spin 
correlator has been derived and compared with theoret¬ 
ical predictions, supporting that the low-energy effective 
theory is the SU( 2 )i WZW model with free boundary 
condition. We also characterize the full spectrum of the 
SU(n) open Haldane-Shastry models and determine the 
twisted Yangian generators responsible for the highly de¬ 
generate multiplets in the energy spectrum. The present 
infinite MPS with open boundaries is readily applica¬ 
ble to any boundary CFT for finding their lattice dis¬ 
cretizations. As an outlook, we expect that the infi¬ 
nite MPS with OBC could be very useful for proposing 
trial wave functions for single-impurity Kondo problems, 
where boundary CFT are known [SS] [39] to play an im¬ 
portant role. 
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Supplemental Material 


Inhomogeneous open Haldane-Shastry models 

In this Section, we provide details on the derivation of the spin-1/2 inhomogeneous open Haldane-Shastry model 
and its SU(n) generalization. 

To construct the spin-1/2 inhomogeneous open Haldane-Shastry model, we use the operators annihilating the 
spin-1/2 open infinite MPS 




(1) 


to build a positive semidefinite operator 

= -^ E + Wi-k){Sj ■ Sk - 2ieabcStS]Sl + SabcSadeS^StS^ISt) 

— g ^ ^ T T ’ ^k ^j^k abc^abe^j ^k ^^abc^ade^bdf 


= -^ E + Wik){Sj ■ Sk - ieabcStS]Sl) 

j,k{j^i) 

= +Wfjf{- + Si-Sj) - - E (Wij+Wij){wik+Wik){Sj ■ Sk), (2) 


where we have used w*j = -Wij, S’^S'^ = \6ab + ^SabcS", Sabc^abd = 2(5cd, and Sabc^adeSbdf = Scef- Then, we obtain 


E(A“)U“ = + Wi-j)^(^ + -5* ■ ^ E ( E + Wij){wik + Wik) I {Sj ■ Sk). 

i,a jj^k J 

The following cyclic identity is the key for simplifying ih 

{Wij + Wij){Wik + Wik) + (Wji + Wfi){Wjk + Wjk) + {Wki + Wki){Wkj + Wkj) = 4. 


By using this identity, we obtain 


(3) 

(4) 


E {Wij + Wij){Wik + Wjfe) = E + Wj^) - (Wki + Wki){Wkj + Wkj)] 

i{^3^k) ii7^j,k) 

= 4(L - 2) - {wjk + Wjk) E E 

K^j,k) ii7^j,k) 


= 4{L - 2) + 2{w]k + wA) - {wjk + Wjk) 


+ E 




-{Wkj+Wkj) 


Wkk + E 

i(jik) 


+ Wjj{Wjk + Wjk) + WkkiWkj + Wkj) 


= {4L - 6) + 2{w% + + Wjk{cj - Ck) + Wjk{cj + Ck), 


(5) 


where we have defined Cj = wjj + and have used Wjj{wjk + + w^j.{wkj + 'a^/cj) = 2 (the latter 

can be easily proved by using the cyclic identity WijWik + WjiWjk + WkiWkj = 1). 
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By substituting <§ into ([^, we arrive at 

Vij + ujij) y- 






^(4L - 6) + + ^1^) + Wjk{cj - Ck) + + Ck) {Sj • 5/^) 




1 1 

+ 


\Zi-Z- 


-.12 


8 E 

“ - \ E(^b + + 2A^ 


Wij{Ci - Cj) + WfJCi + Cj) 


12 


(5i • 5,) 


(6) 




where we have used wfz = 1 — i— 


Then, the spin-1/2 inhomogeneous open Haldane-Shastry model is defined by 


= E 






Wij{Ci - Cj) -\-Wfj{Ci +Cj) 

12 




( 7 ) 


whose ground-state energy E is given by ^ ^ ~ • 

The derivation of the SU(n) inhomogeneous open Haldane-Shastry model follows the similar steps for the spin-1/2 
case. The operators annihilating the SU(n) infinite MPS are given by [Tl|2] 


^ 2(/+^l) E + wqWj + {^^dabc + ifabc)t\^ 


(8) 




where dabc and fabc are the SU(n) totally symmetry tensor and the totally antisymmetric structure constant, respec¬ 
tively. Similar to the spin-1/2 case, we consider the positive semidefinite operator 

E(A“)'^A“ = E + + + 

a '' > jM^i) 

= E + 2(n+ 


= - E + ^ij)^ 


(n - l)(n + 2) (n- l)(n + 2) ^ 


4n 


- E + Wik) 


2(n + l) 


^T2,-»-». Tl n h r 


-dabctitjtl 


_2(n + 1)' ^ ' 2(n+ 1)' 

where we have extensively used the identities listed in the Appendix A in Ref. p. Notice that 
y] (Wij + Wfj){Wik + WiJ.)dabct1t]tl 


( 9 ) 


i^j^k 


= o E + Wij){Wik + Wik) + (Wji + Wfl){Wjk + Wjk) + {Wki + Wkz){Wkj + Wkj)]dabcttt’jtl 


i^j^k 


= o E dabctitjtl 


i^j^k 

4 u 2iin? — 4) 2{in? — l){in? — 4) 

= -dabcT'^T^T^ - -^r“r“ + - tL, 

3 n 6n^ 


( 10 ) 
































where tf. Together with ([^, we obtain 

i,a 


{n — l)(n + 2) {n — l)(n + 2) 


An 


2(n + l) 


(ti ■ tj) 


■ ^ {Wij + Wij){Wik + WiJ,) 
i^j^k 


"+'«>&)+ " - 


— 2(n + 2) 

'■ 

2n 


+ 


2 (n + l)^^^ ' 2 (n + l) 

- Cj) + Wfj{ci + Cj) 




•j \ /^j\ A{ti T 1) 

arr^br^c 2(n + _ (n - l)(n + 2) 


dabctftjtl 

(ti • 


3(n + 1) 

(n - l)(n + 2) 
6n 


dahrT^T^r - 


n + 1 
L(6L + n-2). 


4n 




Then, the SU(n) inhomogeneous open Haldane-Shastry model can be defined as 


^ 2(n + 2) 3(n + l)(n + 2) 


dabcT'^T^T^ + + E 


L 


= E 


1 1 

+ 


n + 1 

Wii(ci-Cj)+M;jAci + Cj) 


\Zi-Zi\‘^ \Zi-Z-j\^ 


^J\ 


4(n + 1) 


{ti -tj), 


( 11 ) 


( 12 ) 


whose ground-state energy E is given by £’ = ~ TEi + n — 2). 


Two-point spin correlation function for the type-I spin-1/2 open Haldane-Shastry model 


In this Section, we derive the exact expression of the two-point spin correlation function for the type-I spin-1/2 
open Haldane-Shastry model. 

As we mentioned in the main text, the two-point spin correlation function Cij = {At\Si • ASj|T)/(T|T) satisfies the 
following linear equations: 


Uj — Un 


-a 


Ui — Ui 4 Hz — 




4 Uj — Un 


'^i,j 


(13) 


where Uj = cosdj. Since is a spin singlet, J2j=i = 0) correlator also satisfies 


E - 4 • 


(14) 


For instance, if one wants to determine the correlators involving the first spin, one could write down the L — 1 linear 
equations (relating Cij, j = 2,..., L) in a matrix form: 


/ 1 

1 

1 

1 



/ 1 \ 

" Wi^—U 2 

U2—U3 

W 2 yW 4 

U2-^Ul 

\ 


l-'12 ] 
Ci3 

Cu 


(M 

U3-U2 

Ui^-U 3 

U3-U4 

Us-^UL 



3 

U3-U1 

U4—U2 

U4—U3 

U\—U4 

U4—UL 



“ ~4 

U4—U1 

1 1 

1 

1 

1 

■J 


\Cil) 


1 1 i 

\ UL—U2 

UL—U3 

UL—U 4 

Ui—UL 



^UL—Ui ' 


(15) 


The correlators involving other spins can be solved in a similar fashion. For the moment, we carry out the derivations 
based on (15) for ease of notation and, in the end, extend the final result to the most general case. 

For the type-I case with Oj = p{j — ^), the following sum identity is very useful: 


E^- 

ji^i) 


cos mOj = 


2{L — m) sin Oi sin mOi — cos Oi cos mOi 
2 sin^ Oj 


(16) 
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where m is an integer and m G [0, 2L]. 


For the l-th row in (|15[), we multiply cosm6>/+i and then sum over all the linear equations. By using (16), we obtain 


L r 

E 

j=2 L 


COS mOi + cos mOj 2{L — m) sin Oj sin mOj — cos Oj cos mOj 
cos 0i — cos Oj 


2 sin^ 0 ^ 


C 


li 


3 2(1/ — m) sin 0i sin mOi — cos 0i cos mOi 


sin^ 0i 


where m G [0, 2L]. For m = 0, this yields 


E 


+ 


5 


cos 0i — cos Oj 2 sin Oj 


Cij = -^ 


3 cos^i 


8 sim 0i 


When multiplying (18) by cosm^i and then subtracting with ( p!7| ), we obtain 


L r 

E 

3=2 L 


COS mOj — cos mOi 2{L — m) sin Oj sin mOj + cos Oj (cos mOi — cos mOj) 


cos 0i — cos Oj 


2 sin^ Oj 


Cij = ^ 


3 {L — m) s'mmOi 


sin^i 


Manipulating three consecutive linear equations [taking m — 1, m, and m + 1 in (19)], we arrive at 


E 


(21/ - 2m + 1) 


cos(m + l)6>j cos{m — l)0j 


sin^ 0^ 


{2L-2m-l)- 


sin^ Oj 


which we have verified to hold for m G [0, 2L]. 

In general, the two-point spin correlator satisfies the following equation: 


E 


(2I/-2m + l) 


cos(m + 1)0j ^ ^ ^ cos(m — 1)0j 


sin^ On 


- {2L-2m-l)- 


sim On 


(17) 

(18) 

(19) 


(cos6>i — cos0j)Cij = 3cosm6>i, (20) 


{cos Oi — cos 0j)Cij = 3 cos mOi, (21) 


where m G [0, 2L]. 

In practice, finding the analytical form of Cjj directly from (21) does not seem to be a simple task. Here we 
adopt an approach used in Ref. [3] to determine the analytical form of Cjj for a few finite-size chains, from which a 
well-educated guess helps to solve (21). 


In the hardcore boson basis, the type-I open Haldane-Shastry ground state is written as 

|^)= ^ ^(:ri,...,Xi/2)5+ ■■■5+^JO), 

Xi<...<Xl/2 


( 22 ) 


where 


L/2 

^'(a;i,..., XL/a) = (-1)^'=!®'sin ff (cos6»j:, - cos6»j:J^. (23) 

1=1 l<i<j<L/2 

Here xi,... denote the positions of the hardcore bosons (up spins). 

The norm of ( [2^ is given by 

(^-l^r) = \^{xu...,Xl/2)\‘' 

Xi<...<Xl/2 


E n (cos6>^.-cos6>^.)^ 

^ xi,...,XLi2l=l l<i<j<L/2 





/I 

cos 

COS^ 

COS^ Ox^ 


cos"^”^ 9x^ 

\ 





0 

1 

2 cos Ox:^ 

3 cos^ Ox^ 

■■■ (L 

- 1) cos^“2 9x^ 





L/2 

1 

cos 0:^2 

cos^ 0x2 

cos^ 0x2 


cos^~'^ 9x^ 



1 

(L/2)! 

E 

Xx,...,Xl/2 

sin^ Oxi det 

1=1 

0 

1 

2 cos 0x2 

3 cos^ 0x2 

■■■ (L 

- 1) COS^“2 


, (24) 




1 


cos2 

cos3 


cos^-1 9x^1^ 






VO 

1 

2 cos Oxj^/2 

3 cos^ 9x^,^ 

... (L- 

- 1) cos^-2 9x^1^ 

) 
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where in the last step we have used the Confluent Alternant identity [3] 


n - VjY = det 


/I 

yi 

2/1 

y \ 

0 

1 

22/1 

32/? 

1 

02 

2/1 

2/2 

0 

1 

2^2 

32/i 

1 

Vm 

yli 

y%[ 

Vo 

1 

2^m 



\ 

(M - l)yf 
{M - l)yf 


2/m ^ 


(M - l)y^-V 

Similarly, the unnormalized transverse spin correlator (for i j) can be expressed as 

(^|5+5-|^) = - 


X det 



i E 



,XlI2-i) 

-1)! 

1 

(M 



y-j 

sin Oi sin Oj 

L/2-1 

11 sin^^a;, 

-1 1=1 


- 1)! cosOi — 

COS Oj ^ 

•J Xi,...,Xl/2- 


/I 

cos Oi 

cos^ Oi 

cos^ Oi 

COS-^“^ di \ 

1 

cos Oj 

cos^ Oj 

cos^ Oj 

COS-^“^ dj 

1 

cos Ox^ 

COS^ Ox^ 

COS^ Ox:^ 

COS^“l 

0 

1 

2 cos Ox^ 

3 cos^ Oxi 

(L-l)cos^-2 0,, 

1 

Vo 

cos6>^^/,, 

1 

-1 

2cos6>^^/,_^ 

cos^ ■ • ■ 

cos^-1 

(L - 1) cos^-2 


(25) 


(26) 


For small L, (24) and (26) can be computed by expanding the determinants (with Laplace’s formula). After the 
expansion, the discrete sums over the coordinates can be carried out by using the following identities: 


1 1 

sin^ Ox cos^’^ Ox =- 5 - 

^ ^ r 22^- 

X = 1 


2r 

22r+l \ ^ 


(27) 


and 


sin^ Ox cos^’^^^ Ox = 0, 


X = 1 


which are valid for the type-I case and r = 0, 
Following this procedure, we obtain for L = 4 


^-1 

. , 2 -L- 


(«'|5+5y^') (-1)*-^' sin^isin^,- 


(^1^) 


L cos Oj — cos 0, 


-[2{cosOi — cos6>j) + - {cos 30i — cos30j) 


— - (cos 20i cos 30j — cos 30i cos 20j)]. 
0 


For 1/ = 6, we obtain 

|4^) (—1)*“-^ sin6>^sin6>j 


(28) 


(29) 


6 14 

^ ^ ^ [2(cos^i — cos6>A +-(cos36>i — cos3^A +--(cos56>i — cos5^A 

L cosOi-cosOj^ ^ 15^ 

4 28 

— - (cos 20i cos 30j — cos 30i cos 20j) — ^ (cos 20i cos bOj — cos bOi cos 20j) 


— - (cos 4:0i cos bOj — cos bOi cos 40j)]. 


(30) 
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For L = 8 , we obtain 

1 ^) (^—ly-j sinOisinO^ 


6 14 

, ^ ^ , - — ----- ^\2(cos Oi — cos Oj)-\—(cosSOi — cosSOj)-{ -(cos 56^^ — cos 56>.) 

{^\^) L cosOi-cosOj^ ^ 15 ^ 

154 4 28 

+ -^ (cos 70i — cos 70j) — - (cos 20i cos 30j — cos 30i cos 20j) — (cos 20i cos bOj — cos bOi cos 20j) 

ic/O O 4t) 

308 4 

-(cos 20 i cos 70^ — cos 70i cos 20 n) -(cos AOi cos bOn — cos bOi cos AO^) 

585 V ^ 3 ^ 33 ^ 3 t 3J 


44 4 

— (cos 4:0i cos 70j — cos 70i cos40j ) — — (cos 60i cos 70j — cos 70i cos bOj)]. 

Since |4^) = and |4^) is a spin singlet, we have 


Cij = 


^3{^\StSj 1 ^-) 


(31) 


(32) 


For larger L, a direct computation of (24) and (26) becomes quickly involved. However, from the finite-size results 
(29)-(31), there is an indication that, for general L, the analytical form of the two-point spin correlator Cij is given 
by 


L12 p—1 


Cij = 


3(—1)* sin sin 
L{cosdi-cosej) 


EE gpq[cos{2p — l)0i cos 2q0j — cos 2q0i cos( 2 p — 1)0j]^ 


(33) 


where Qpq has no L dependence and its initial values are readily available from (31). 

By substituting (33) into (21), the well-educated guess (33) indeed solves the linear equation and the general 
expression for Qpq is found to be 


9pq — 


1 

p-l 4m-l 


llm=l 


p = l,g = 0 
p > l,g = 0 . 


4m+l 

o rr^-l 4m-l Tjq 4n-3 ^ > 1 o > 0 
^llm=l 4m+l iln=l An-1 P ^ ^ ^ 


(34) 


Twisted Yangian generators for the SU(n) open Haldane- S hastry model 


In this Section, we provide details on the derivation of the twisted Yangian generators for the SU(n) open Haldane- 
Shastry model. 

For the SU( 2 ) open Haldane-Shastry model, such formalism has already been developed in Ref. [4]. Althought its 
SU(n) generalization is rather straightforward, we present the derivation below for the purpose of being self-contained. 
Following Ref. [4], we introduce an unprojected Hamiltonian 


^ = -E 


ZiZj 


{Zi -ZjY 


{Kij - 1 ) + 


ZiZ. 


-1 


(Zi-Zi ) 




L r 


-E 


hi 


Zi 


{zi - 1)^ 


+ 21)2 


{zi-zpr 


{Ki - 1), (35) 


where the coordinates Zi are viewed as dynamical variables, the coordinate permutation operators Kij, Kij^ and Ki, 
when acting on the coordinates, yield KijZi = zjKij^ RijZi = zJ^Kij^ and KiZi = zy^Ki^ and the constants hi and 
62 will be specified below. 

We also define a projection operation tt which replaces the operators Kij and Kij by the SU(n) spin permutation 
operator Pij = 2ti • and Ki by the identity operator once they have been moved to the right of an expression. 

In the simplest case with only one of these operators, we have 

n^Kij) = 7i{Kij) = Pij, (36) 

^{Ki) = 1 . (37) 

If there are multiply coordinate permutation operators Kij and Kij present, the rule of the projection operation is to 
insert a designed product of SU(n) spin permutation operators (which itself should be an identity, e.g., PikPijPijPik = 
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1) into the expression and then replace each combined product PijKij (appearing to the right of an expression) by 
an identity, e.g., 

7T{KijKik) = 7T{PikPijPijKijPikKik) = PikPij‘ (38) 

After the projection operation, the coordinates are not dynamical any more. Then, the projected Hamiltonian is a 
pure SU(n) spin model 


H = 7r{H) 

= -E 

= E 


{Zi -Zjf- 

1 


{Pij - 1) + 


ZjZ 


-1 




+ 




(zi-zpr 

(Pij - !)• 




(39) 


In Ref. [4], it has been shown that the projected Hamiltonian is integrable, if the lattice coordinates correspond to 
the three uniform cases (see Fig. 1 in the main text) and the constants bi and 62 in (35) are given by (i) type-I: bi = 0 
and 62 = 1; (ii) type-H: 61 = 0 and 62 = 3; (hi) type-HI: 61 = 62 = 1. Notice that the three projected Hamiltonians 
( [^ , after subtracting a constant, just correspond to the open SU(n) Haldane-Shastry model [Eq. (7) in the main 
text]. 

The integrability becomes manifest by introducing the Dunkl operators 


“ E E .. 


-K, 


j{>i) 




E 


i(<d 




bi- 


Zi Z- 


3T K- 


(40) 


which are mutually commuting, [di^dj] = 0 Vi, j, and all commute with the unprojected Hamiltonian, [di^H] = 0 Vi. 
After introducing an extra n-dimensional auxiliary Hilbert space (denoted by “0”), the SU(n) monodromy matrix 
T{u) can be defined as 


T{u) = TT 


no 


iO 


u — dj 


1 


+ 62 1 


u 


no 


iO 


u di 


(41) 


u=l ^ ^ ^ i=L 

which is a n X n operator-valued matrix function of the spectral parameter u. Actually, it is a generating function of 
conserved charges, [T{u),H] =0 \/u. By using the Taylor expansion ^ ^ and implementing 

the projection, one obtains formally the following expression: 


n^-l 


n ^-1 


T{u) — 1 -\— ^ ^ tg (8) Tf j H —2 [ ^0 ^ ^2 V ^ ^ ^ '^2 ] V • • • , 


(42) 


a=l 


a=l 


where and (a = 1,..., — 1 and /i = 1,..., 00) are conserved charges for the SU(n) open Haldane-Shastry 

model, = 0. For the monodromy matrix (41), the conserved charges in the first- and secord-order 

expansions in 1/u are trivial (such as T^, dabcT^T^i ^ etc). In the third-order expansion, we obtain, after a 

tedious but straightforward calculation, the following nontrivial conserved charge: 


= E+T'lw^fco) - 72 E + '^jk)iwij - wfj)tiPjkPij 

k 


(43) 


where 71 and 72, for the three uniform cases, are given by (i) type-I: 71 = 0, 72 = (ii) type-H: 71 = 0, 72 = 
(hi) type-HI: 71 = 1, 72 = respectively. 
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